Gravitational dynamics in s + 1 + 1 dimensions II. Hamiltonian theory 
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We develop a Hamiltonian formalism of brane-world gravity, which singles out two preferred, 
mutually orthogonal directions. One is a unit twist-free field of spatial vectors with integral lines 
intersecting perpendicularly the brane. The other is a temporal vector field with respect to which we 
perform the Arnowitt-Deser-Misner decomposition of the Einstein-Hilbert Lagrangian. The gravi- 
tational variables arise from the projections of the spatial metric and their canonically conjugated 
momenta as tensorial, vectorial and scalar quantities defined on the family of hypersurfaces contain- 
ing the brane. They represent the gravitons, a gravi-photon and a gravi-scalar, respectively. From 
the action we derive the canonical evolution equations and the constraints for these gravitational 
degrees of freedom both on the brane and outside it. By integrating across the brane, the dynamics 
also generates the tensorial and scalar projection of the Lanczos equation. The vectorial projection 
of the Lanczos equation arises in a similar way from the diffeomorphism constraint. Both the gravi- 
ton and the gravi-scalar are continuous across the brane, however the momentum of the gravi-vector 
has a jump, related to the energy transport (heat flow) on the brane. 



I. INTRODUCTION 

The Hamiltonian theory of general relativity is based 
on the Arnowitt-Deser-Misner (ADM) decomposition [l[ , 
i.e. slicings of the four dimensional (4d) space-time man- 
ifold with a family of spacelike three-dimensional (3d) 
hypersurfaces. 

The slicing method with respect to a timelike hypersur- 
face (the brane) has been also proven efficient in brane 
world models which interpret our traditional 4d space- 
time as a brane embedded in a 5-dimensional (5d) man- 
ifold, the bulk (for a review see 0)- In the 4+1 decom- 
position of the bulk, various geometrical projections of 
the 5d gravity on the brane were used together with the 
Lanczos-Sen-Darmois-Israel (LSDI) matching conditions 
in order to derive the effective Einstein equation 
on the brane 0]-@|. In these approaches the 4d space- 
time manifold was treated covariantly. In the covariant 
approach however it is not obvious how to describe time 
evolution, and introduce the very concept of the (many- 
fingered) time function, which is the first step in devel- 
oping a Hamiltonian formulation of the brane world sce- 
nario. 

In a previous paper 0] (hereafter quoted as Paper I), 
we have discussed a two-fold slicing of the bulk, both 
with respect to a timelike and a spacelike foliation. One 
of them was meant to contain the brane, the other the 
constant time slices St (referring to £=const) of the bulk, 
necessary for Hamiltonian evolution. The result was a 
(3 + 1 + l)-decomposition of the bulk, which resulted in 
a break-up of the brane into space and time. The formal- 
ism was given in fact for the more generic case of s + 1 + 1 
dimensions, therefore it is applicable for general relativis- 
tic situations as well, whenever any space-like direction 
should be singled out for one reason or another. 

In this framework the gravitational degrees of freedom 
are represented by new variables. These are the spatial 
metric g a b on the brane, the shift vector M a and lapse 



function M associated with the brane normal; together 
with the tensorial, vectorial and scalar projections of the 
extrinsic curvature associated with the temporal normal. 
These can be thought of as gravi-tensorial, gravi- vectorial 
and gravi-scalar degrees of freedom, the extrinsic curva- 
tures representing their generalized velocities. 

In Paper I we have given the decomposition of all ge- 
ometric quantities with respect to the two normals, to- 
gether with the evolution equations of the gravi-tensorial, 
gravi- vectorial and gravi-scalar degrees of freedom in the 
velocity phase-space. A set of (Lagrangian-type) con- 
straint equations on these variables were also given. The 
latter should be imposed on any set of initial data on a 
spacelike section of the brane. 

However Paper I did not contain any reference to 
canonical momenta and the phase-space of gravitational 
degrees of freedom. Neither did it address the variational 
principle leading to the Hamiltonian evolution equa- 
tions, together with the super-Hamiltonian and super- 
momentum (diffeomorphism) constraints of brane-world 
gravity. 

It is the purpose of this paper to discuss these topics, 
together with the regularization issues raised by the pres- 
ence of the distributional matter sources on the brane. 
In Sec. II we give a brief, somewhat technical account of 
those results presented in Paper I which are needed for 
our present considerations. 

Sec. Ill contains the variational principle for vac- 
uum gravity, with concomitant ADM and brane-bulk 
decompositions. We define here the momenta associ- 
ated with the gravi-tensorial, gravi-vectorial and gravi- 
scalar degrees of freedom and perform the Legendre- 
transformation. Extrcmizing the action with respect to 
the gravitational phase-space variables give the Hamilto- 
nian evolution equations in terms of variables defined on 
the t =const space sections of the brane. Extremizing 
with respect to the lapse iV a and shift N (both related 
to the temporal evolution) gives the super-Hamiltonian 
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and super-momentum constraints of brane-world gravity. 
Both the shift N a and super-momentum TL^ have only 
s components, as the off-brane component of the shift is 
suppressed in order to obey Frobenius' theorem, i.e. to 
let the brane exist as a hypersurface. 1 

Matter sources both in the bulk and on the brane 
are discussed in Sec. IV. We show here the following 
commutativity property: the (s + 1 + l)-dccomposition 
of the energy-momentum tensor obtained by a covariant 
variation leads to the same result as the non-covariant 
variation with respect to the gravitational variables aris- 
ing from the (s + 1 + l)-dccomposition. We also give 
here the Hamiltonian evolution in the presence of mat- 
ter. Similarly as the evolution equations in the velocity 
phase-space given in Paper I, these Hamiltonian evolu- 
tion equations contain distributional type sources. 

We remedy this situation in Sec. V by a regulariza- 
tion procedure across the brane. We consider an interval 
(— XO; Xo) across the brane (lying at x = 0) and integrate 
the supcr-Hamiltonian and super-momentum constraints 
of the total system (gravity + matter) . We also integrate 
the Hamiltonian evolution equation of the gravi-tensorial 
momentum. Finally we take the thin brane limit xo - *• 0. 
This method gives the scalar, vectorial and tensorial pro- 
jections (with respect to S t ) of the Lanczos equation. 

Sec. VI contains the concluding remarks. We include 
three Appendices, the first containing the (s + 1 + 1)- 
decompositions of various quantities needed throughout 
the paper. In the second Appendix we schetch the deriva- 
tion of the projections of the Lanczos equation in the ve- 
locity phase-space. The third Appendix proves another 
commutativity property. While in the main paper we 
obtain the results by an s + 1 + 1 ADM decomposition 
prior to applying the variational principle; the same re- 
sults can be regained by first performing an (s + 1) + 1 
ADM reduction, then extremizing the action and finally 
by a further(s + l)-split. 

Notation. — All quantities defined on the full (s + 2)- 
dimensional space-time and on the (s + l)-dimensional 
brane carry a distinguishing tilde and the prefix (s + 1), 
respectively. A hat distinguishes quantities defined on 
the spatial (s + l)-leaves. Quantities defined on the spa- 
tial sections of the brane carry no distinctive mark. For 
example, the metric 2-forms are denoted g, g and g, 
respectively, while the corresponding metric-compatible 
connections are V, D and D. Latin indices represent ab- 
stract indices running from to (s + 1). Vector fields 
in Lie-derivatives are represented by boldface characters. 
For example £\-T denotes the Lie derivative along the in- 



While in the s + 1 ADM decomposition the congruence d/dt gen- 
crating time-evolution should be twist-free in order St to exist, 
there is no such condition, that an (s — l)-parameter subset of 
the congruence should form a temporal hypersurface. It is ex- 
actly this condition, which should be obeyed in brane-worlds, 
leading in the simplest case to the suppression of the (s + 1)* 
component of the shift. 



tegral lines of the vector field V a . A dot denotes a deriva- 
tive with respect to t, such that T = (d/dt) T = £ d / dt T. 



II. THE (s + 1 + ^-DECOMPOSITION OF 
SPACE-TIME 

In order to develop the Hamiltonian theory of brane- 
world gravity from the viewpoint of an observer on the 
brane, a double decomposition of the bulk is necessary. 
We have to separate both the normal to the brane l a 
and the temporal normal n a , with respect to which the 
ADM decomposition will be carried out. These vector 
fields obey n a n a = —1, l a l a = 1 and n a l a = 0. We 
introduce a coordinate x transverse to the brane and we 
foliate the (s + 1 + l)-dimensional space-time by (s + 1)- 
dimcnsional leaves of constant x> the brane being at x = 
0. The off-brane evolution is defined along the integral 
lines of 

(£) = Ma+Mia > w 

where M a and M represent the shift vector and lapse 
function associated with the off-brane evolution. They 
characterize the off-brane sector of gravity. 

The temporal evolution is along the integral lines of 

(I)" = Nn a + N\ (2) 

where N a and N are the familiar shift vector and 
lapse function related to the foliation of the (s + 1 + 1)- 
dimensional space-time by the (s + l)-dimensional spa- 
tial leaves t =const. The coordinate t represents time. 
Note that neither (d/dx) a has a component along n a , 
nor (d/dt) a along l a . The first condition guarantees that 
(d/dx) a generates a purely spatial displacement. The 
second one is the simplest gauge-choice, which ensures 
that the Frobenius theorem holds, thus l a is hypersurface- 
forming (see Appendix C of Paper I) . 

The metric g ab of the (s + 1 + l)-dimensional space- 
time, decomposed with respect to l a and n a is 

g a b = g a b - n a n b + l a l b , (3) 

where g a b is the induced metric on the s-dimensional 
leaves S tx , which represent the intersection of the 
X=const and i=const leaves. The components of (s + 2)- 
dimensional metric in the coordinate system adapted to 
the coordinates (t,x,x a ) can be written as 

ds 2 = -(N a N a - N 2 )dt 2 + 2N a dtdx a + N a M a dtd X 
+g ab dx a dx b + 2M a dx a d X + (M a M a + M 2 )d X 2 . (4) 

Then, as discussed in Paper I, the gravitational sector 
is described by {g ab , M a , M, N a , N} obeying g ab n a = 
g ab l a = M a n a = M a l a = N a n a - N a l a = 0. The 
above set has one variable less then the number of vari- 
ables contained in g ab . This is because time-evolution is 
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restricted to proceed along x=const hypersurfaces. (Ac- 
cording to Eq. ([21) there is no off-brane component of the 
shift.) 

In Paper I we have introduced two types of extrinsic 
curvatures, related to the normal vector fields n a and l a . 
Each of them could be further decomposed with respect 
to the other normal into tensorial, vectorial and scalar 
projections. These give two kinds of second fundamental 
forms of the leaves T, tx 

K ab = g c a 9 d b^cn d , L ab = g c a g d b^cld , (5) 

two normal fundamental forms 

JC a = g b J c V c n b , C a = ~g b a n c V c l b (6) 

and two normal fundamental scalars 

AC = l a l b V a n b , C = n a n b V a l b . (7) 

From the symmetry property of the extrinsic curvatures 
and n a l a — the condition: JC a — C a follows. 

It was shown in Paper I that L ab and C can be 
expressed in terms of ^-derivatives and the covariant 
derivatives D a associated with g ab as 



L ab = ^-A^t-W^.M, 



2M \ dx 



C = -J-(™-M"D a N 
MN V d X 



(8a) 
(8b) 



Thus, L ab and C are related to the spatial derivatives of 
the gravitational variables. By contrast, the quantities 
K ab , )C a and K. give the time evolution of g ab , M a and 
M, respectively: 



K„ 



K a = 



1 / dg ab 



2N \ dt 



2D (a N b) 



1 (dM a dN a 



2MN\ dt 



dx 



-A fD b N a -N b D b AI a ) (9b) 



MN V dt 



(9c) 



therefore they are velocity- type variables. 

Thus gravitational dynamics in the velocity phase- 
space can be given in terms of {g a b, M a , M; K abl K. a , K.}. 
Time-evolution of {g ab ,M a ,M} is expressed by Eqs. 
(|9]), while equations representing the time-evolution of 
{K ab ,K, a ,/C} were given by Eqs. (67) of Paper I . 

In the next section we present the Hamiltonian formu- 
lation of brane-world gravity, by passing to the momen- 
tum phase-space. 



III. HAMILTONIAN EVOLUTION OF VACUUM 
GRAVITY IN THE BULK 

Vacuum geometrodynamics (without cosmological 
constant) arises from the Einstein-Hilbert action 



S G [g ab ] = J d s+2 xC G = J d s+2 x^R 



(10) 



The (s + 1 + l)-decomposition of the (s + 2)-dimensional 
scalar curvature and of the metric determinant were de- 
rived in Paper I as Eqs. (61) and (B2), respectively. The 
vacuum gravitational Lagrangian density becomes 

C G [g ab ,M a ,M,K ab ,JC a ,JC-N\N] = 
NMy/g(R - L ab L ab + L 2 - 2CL) + 2^§D a ND a M 
+NM^(K ab K ab -K 2 + 2K. a K. a - 2KK) 
-2V a [NM y /g~(a a - X a - Kn a + Ll a )] . (11) 

Here a a = n b V b n a = N- 1 D a N~Cl a and A a = l b V b l a = 
—M~ 1 D a M + K.n b are acceleration-type quantities (the 
curvatures of the normal congruences n a and l a ). By 
transforming (the double of) the terms quadratic in the 
extrinsic curvatures in the Lagrangian density (jlip into 
time-derivative terms cf. Eqs. ([9]) we obtain the following 
advantageous expression: 

C G [g ab ,M a ,M,K ab ,JC a ,lC-N a ,N] = 
JgM [K ab — {K + K) g ab ] g ab + 2J~gK a M a 
-2^/gKM - NHf - N a H G 

-2y/gD a [M{D a N + N b K ab ) + NCM a - N b JC b M a ] 
+2— [y/g (NC- N a K, a )] + 2- [M*Jg(K + K)] , (12) 

where we have denoted 



Hi = 



n, 



~^g[M(R + L 2 -3L ab L ab ) 
~2g ab (d/d X - £m)L ab - 2D a D a M 
+M{K 2 - K ab K ab - 2tC a tC a + 2KK)] , (13a) 
-Jg{D b [MK\ - Mg b a (K + AC)] + KD a M 



+MK a L + {d/d X - £m)K a } 



(13b) 



(Dots represent the time-derivatives d/dt). The advan- 
tage of writing the Lagrangian density in the form (|12|) 
is that it contains explicitly the Liuville form, with the 
right coefficient to drop out when performing later on 
the Legendre transformation. The cofactors TC G and H G 
of the Lagrange multipliers N and N a are the super- 
Hamiltonian and the super-momentum constraints of 
vacuum gravity, also obtainable as projections of the Ein- 
stein tensor, cf. Appendix A: 



H G 



-2M^n a n b G ab 



H G = ^2My^g b a n c G, 



be 



(14a) 
(14b) 



Since the (s + l)-th component of the shift was set to 
zero, the super-momentum contains only the components 
corresponding to brane spatial diffeomorphisms. 

Now we define the phase space of the brane vacuum 
gravity as the set of canonical coordinates and canoni- 
cally conjugated momenta, 



{g A ;TT A } := {g ab ,M a ,M;TT ab ,p a ,p} , 



(15a) 



by introducing the notation g A = {g ab ,M a ,M} and 
tt a = {ir ab ,p a ,p} with the multi-index A = 1,2,3 
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as a condensed notation for the gravi-tensorial, gravi- 
vectorial and gravi-scalar degrees of freedom. The mo- 
menta canonically conjugated to the field variables g a b, 
M a and M are 



7T 
Pa 
P 



dC G 
dtjab 
dC G 

dM a 
d£ G 



M^g[K ab -{K + K)g ah ] , (16a) 
= 2^/gK a , (16b) 



dM 



v = -2^/gK . 



(16c) 



Inverting Eqs. (|T6|) with respect to the second and 
normal fundamental forms and the normal fundamental 
scalar, we obtain 



K ab = 

K. a = 
K = 



1 



1 



p" 



P 



2^9 
s-1 



-P 



2s^/g r sM^g 



g a b , (17a) 
(17b) 
(17c) 



We insert these formulae into the Lagrangian (jTTj) . 
so that the Einstein-Hubert action takes the "already 
Hamiltonian form" : 



S G [g A ,n A ;N a ,N} = 



dt 



dx 



' —oo J— OO *J Z->tx 

x(7r A g A - NHfiQA,^] - N a n G [g A ^ A ]) 
—ST t= ±oo — STgs tx + ST x= ± oa . (18) 

The extremization of the action (fl8|) with respect to the 
lapse and shift gives the constraint equations 



SS G 
' SN 



n, 



ss G 

'SN a 



(19) 



In detail, the super-Hamiltonian constraint Ti. G and the 
super-momentum constraint Ti G of vacuum gravity, writ- 
ten in terms of the canonical data, are: 



-^g-[M{R-L 2 + 3L ab L ab ) 

-2g ab (d/d X - £w)L ab - 2D a D a M] 

pir 



1 L ^ 

Mjg \ s 



s- 1 



s^/g 



M f„a 



P Pa + 



2s 



H G [g A ,n A ] = -2D b Ti b a - (d/d X - £m)Pa 
+pD a M . 



(20a) 



(20b) 



The contributions ST t= ±oo, STgs tx and ST x= ± oc in 
the action (fl~8|) denote a collection of surface terms on 
t — > ±oo, dYi tx and x ~ > ±oo , respectively. Their con- 
tribution can be compensated by adding surface terms to 
the action. The contributions dY* tx and x ~ > ioo do not 



have analogues in the standard ADM decomposition; and 
they come from partial integrations meant to transform 
the set of variables (N, N a ) into Lagrange- multipliers in 
the action. For completeness, we enlist these terms: 



ST, 



t — ±00 



OO f oo 



dt d X d s x 

— oo J — OO -/Sty 



M 



(21a) 



ST d 



oo />oo /• 

dt d X d s xD a [N b Tr ab 

— oo J — oo **^>tx 

+ V5 (NCM a + MD a N)] , (21b) 

oo poo 



ST x=±oc = 2 dt d X d s x 



oo J — oo 



dx 



- N 



a Pa 



(21c) 



The Hamiltonian of vacuum gravity can be written 
as the smearing of the super-Hamiltonian and super- 
momentum constraints with the lapse function and shift 
vector, 



H G [N] = H^[N] + H G [N a ] , 



(22) 



where 



H G [N] = fd X f dx s N(x,x)H^(x,x), (23a) 

H G [N a ] = fdxf dx s N a (x,x)H G (x, X ) .(23b) 

The Poisson bracket of any two functions fix, x; g A , tt a ] 
and g(x, x; g Al ir A ] on the phase space is defined as 

{f(x,x),h(x', X J )} = 

„ 6f(x, X ) Sh(x', X ') 



d X " / d s x 



Sg A (x",x") 5ir A (x",x") 
, // f , s // Sf{x,x) Sh(x',x') 

X L x X 5v A (x",x")5g A {x",x") ' 

(24a) 

which provides the following non- vanishing commutation 
relations for the Poisson brackets of the canonical vari- 
ables: 

{g ab (x,x)^ cd (x',x')} = S c (a S d b) 6(x,x;x', X ')(?te) 
{M a (x, X ),Pb(x',x')} = S a b S(x,x;x',x') , (25b) 
{M(x,x),p(x',x')} = S(x,x;x',x) ■ (25c) 
or with the condensed notation 

{gA(x,x),K B (x',x')} = 3 B aS(x,x;x',x) ■ (26) 

Here S B A = {S c ia S d b) , S\, 1}, S(x,x;x', X ') - 
5(x, x')S(Xi x') an d the Dirac delta distribution behaves 
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under coordinate transformations as a scalar in its non- 
primed arguments, and as a scalar density of weight one 
in its primed ones. 

The constraints represent restrictions on the initial 
data and their Poisson brackets "close" according to the 
Dirac algebra: 

{nlix, x),nf(x', x ')} = g ab (x, x )n G (x, X )S, b (x, x; x', X ') 

- (x X <-> x'x') , (27a) 

{nf(x, x),n°(x', x ')} = nf(x, x )s, a (x, X ; x', x ') 

+n G a (x, X )S(x,x;x',x') , (27b) 
{H G (x, X ),H G (x', x')} = H G 8, a {x, x; x', X ') 

— (axx <-> bx'x') ■ (27c) 

Time evolution of the canonical data is generated by 
the Hamiltonian of the system as 



cab/ A A\ 
6 (7T ,7T ) 



1 



a be „^ab 
7T r 7T — — 7T7T 



1 



„cd „2 \ „ab 



M „ Qb f 1 



g A (x, X ) = {g A (x,x),H^[N}} = 



TT A (X, X ) = {7T A (X, X ),H G [N}}-. 



G[ATU _ SH G [N] 



5ir A (x,x) 
5H G [N] 
$9a(x,x) 



, (28a) 



.(28b) 



Computation gives the dynamical equations for the gravi- 
tensor, gravi-vector and gravi-scalar degrees of freedom: 



9ab 

M a 
M 



2N 

+£wgab , 

MN „ dN a 
^P a + -5— + £nM 

MN ( s - 1 



(29a) 
(29b) 



sV _ v 2 p--«) +£ „M (29c) 



1 c 

"xp - t;PcP 



s - 1 



2^g \sM 2 

1 ( X ab , ^La h 



4s 



V s 



S(^,^) 



1 ^ab 1_2 
TTafcT 7T 



s - 1 

- . + "1 — 3 



(31a) 



(31b) 



while V ab , V a and V represent the tensorial, vectorial 
and scalar projections of the force term of the (s + 1)- 
dimensional scalar curvature potential: 

V ab (g A ) = -My/g (G ab + 2L ac L b c - LL ab ) 

-^M^(3L cd L cd ~L 2 )g ab 

+V9 {g ac 9 M - 9 ah 9 cd ) (d/d X - £n) L cd 
+y/~g~(D a D b M - g ab D c D c M) , (32a) 
V a (g A ) = -2y/g-(D b L b a -D a L) , (32b) 
V(g A ) = Vg{R + L ab L ab - L 2 ) . (32c) 



The evolution equations (|30|) together with the con- 
straints (|20p are equivalent with the vacuum Einstein 
equations in the bulk. Once the constraints are obeyed 
at some instant of time, the dynamical equations assure 
that they will continue to be satisfied later on. 



and 



■ab 



NS ab + NV ab - NMy/g~£(L ab - Lg ab ) 
+yfg~{D a D b N - g ab D c D c N - g ab D c ND c M) 
+Vgg ab (d/d X ~ £ M )(N£) + £ N ir ab , (30a) 
Pa = NV a -2^[L b a D b N + D a (NC)] 

+£^p a , (30b) 
p = NS + NV - 2yfg(L£ + D a D a N) 

+£np ■ (30c) 

Here S ab and S denote the tensorial and scalar projec- 
tions 2 of the geodesic spray of the DeWitt super-metric 
0: 



2 Note that no vectorial projection S a is present in p a . However 
as time-derivatives and index raising do not commute, there will 
be vectorial kinetic-type terms in p a . 



IV. HAMILTONIAN DYNAMICS WITH 
MATTER SOURCES 

A. General considerations 

The basic scheme of the Hamiltonian formulation does 
not change if we couple matter fields to gravity. We only 
have to enlarge the phase space with the canonical vari- 
ables of the matter sources. The total action describing 
the system is 



S = S G [g ab ]+2n 2 S M [g ab ^ 



(33) 



where k 2 is the gravitational coupling constant in s+2 di- 
mensions and we assume the matter action S M [g ab , ^i] 
contains the metric only in non-derivative terms. This 
assumption is obeyed for all physically relevant matter 
fields and it assures that (a) the vacuum gravitational 
momenta (|16[) remain unchanged in the presence of mat- 
ter; (b) by performing the Legendre transformation, the 
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matter contribution to the total Hamiltonian is just mi- 
nus its contribution to the total Lagrangian. 

Extremizing the matter action with respect to the mat- 
ter fields yields their evolution equations. Extremiz- 
ing with respect to the metric yields, by definition, the 
energy-momentum tensor 



rpab 



2 SS 



M 



f) 



(34) 



We can perform the (s + 1 + l)-decomposition of the 
matter Lagrangian density, without specifying it explic- 
itly. For this we write the energy-momentum tensor as 



fab = (f cd g a c g b d ) + (f cd n c n d j n~n" 
+ (T cd l c l^ l a l b - 2 (T cd gi a n d ) 
+ 2 (f cd g^ a l d ^ l b) - 2 (f cd n c l d ^ n {a l b) , (35) 

and we note that the decomposition Eq. ^ of the metric 
gives 



Sg a b = 8g ab - 2ri( a 5n b ) + 2l( a 6l b ) 



(36) 



The variation of the matter action with respect to the 
metric (after partial integrations) results in 



5~ g S 



M 



d S+2 



ss M 

X — OQnb 



Sg a b 



d s+2 xNM^g^(f cd g a c g b d ) 6g ab 

f cd g> d ~ (f cd n c n d ^ n a + (f cd n c l d ^ Z a ] Sn a 
f cd g a c l d - (f cd nj d ) n a + (r cd l c l d ) l a ] Sl a 
+ (f cd g bc n d ) Sn b ~ (f cd g bc l d ) 5l b ] . (37) 



We would like to replace the variations 5n a , Sl a , Sn a , Sl a 
by the variation of our chosen gravitational variables. As 
d/dt and djd\ are directions unaffected by the variation 
of the metric, from Eqs. ([T]) and ^ we obtain 



Sn a 

sr 



SN 
— ? 

N 



8N a 



fSM in SM a 

1 H 

\ M M 



(38a) 
(38b) 



As, cf. to Paper I , the dual bases are related as {dt) a = 
n a /N and (dx) a = l a /M, the variation of the co-vectors 
also arises: 



Sn a 



SN 
AT' 



Sla 



SM 



(39) 



and we obtain the desired formula: 
SgS M [g ab ,M a ,M;N a ,N; *] = 

Jdtjdx ^ <Px^g{™ (f cd g a c g b ) Sg ab 

+N (f cd l c l d ^ SM + N (f cd g b l d ^ g ab SM a 

-M (T cd n c n d ) SN - M (f cd g b n d ) 3ab <5A Q }(40) 

(Due to the non-derivative coupling of matter to gravity 
there are no momenta dependencies.) The result ([4*0]) 
shows that extremizing the total action with respect to 
the lapse function N and the shift vector N a (similarly 
to the prescription (|19p ) gives the super-Hamiltonian and 
super-momentum contribution of the matter fields 

(41a) 
(41b) 

The super-Hamiltonian and super-momentum con- 
straints of the total system can be written as 

n x = nf + nf m o , 



Uf = 2K 2 M^/g(T cd n c n d ) 
Hi 1 = 2K 2 M^(f cd g b n d )g ab . 



= H a +H a 



M 







(42a) 
(42b) 



with the vacuum and matter contributions given by Eqs. 
([2H)l and Eqs. (|4"Tj) . Here w denotes weak equality (hold- 
ing on the constraint surface in the phase space). 

In what follows, we discuss the canonical equations in 
the presence of matter. Due to the non-derivative cou- 
pling, Eqs. (f2U)l remain valid in the presence of matter 
(since 5S M /Sn A = 0). However the evolution of the mo- 
menta receive additional contributions. Due to remark 
(b), the matter contributions to the left hand side of Eqs. 
(|30p can be found by extremizing the action with respect 
to the dynamical variables g ab , M a and M: 

SS " ^(f'VM. (43a) 



Sg a b 
SS M 
5M a 
SS M 
SM 



= A\ 7/ ( /"',/,„/ 
= N^g(T cd lJ d 



(43b) 
(43c) 



The dynamical equations for ir A with the contributions 
flU} take the form 

7r Qb = NS ab + NV ab - NM^C{L ab - Lg ab ) 

+^g(D a D b N - g ab D c D c N - g ab D c ND c M) 
+V99 ab (d/d X - £m)(NC) + £ N ir ab 
+n 2 NM^(f cd g a c g b d ) , (44a) 
Pa = NV a NV a - 2^g[L b a D b N + D a (NC)] 

+£ NPa + 2£ 2 AVs (f bc g ab lc) , (44b) 
p = NS + NV — 2^/g(L£ + D a D a N) 

+ £KP + 2K 2 N^g(f ab l a l b ) . (44c) 
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These formulae are valid for any matter source cou- 
pled non-derivatively to the (s + 2)-geometry. The mat- 
ter contributions can be further specified, once the 
energy-momentum tensor (or equivalently, the matter 
Lagrangian) is known. 



B. Brane- world scenario 

In the brane-world scenarios the stress-energy tensor 
is decomposed as 

T ab = Uab + [-X (s+1) g a b + {s+1) T a b}S(x) , (45) 

where the regular part H a b represents the non-standard 
model bulk sources, while the distributional term con- 
tains the brane tension A and the energy-momentum ten- 
sor of standard model matter field localized on the brane. 
First we decompose the bulk energy momentum II a b with 
respect to the off-brane normal l a as: 



1 



IU = {s+1) V ab + 2l (a (s+1) ¥b) + lahV • (46) 



Then we decompose ^ s+1 ^V a b and ( s+1 )*p o further with 
respect to the time-like normal n a as: 



(s+l)*P o 



V ab + 2n (a V b) + n a n b V, (47a) 
Pa + Pn a , (47b) 



where V a bn b = V a n a = = P a ?i Q . (Note that g ab V a b ^ 
V, unless ^ s+1 ^Vab happens to be traceless.) 

The brane contribution can be algebraically decom- 
posed as 

(s+1) T Qb = pn a n b + Pcjab + IU + 2n (a Q b) (48) 

with respect to the 4- velocity n a of the fluid, and in terms 
of the energy density p, isotropic pressure P, anisotropic 
stresses U a b and the energy transport (heat flow) Q a 
(here g ab U ab = n a U ab = n a Q a = 0). 

Choosing the 4-velocity of the fluid as normal to the 
spatial slices does not restrict the arbitrareness of the 
foliation. Indeed, the foliation is given by the form 
(dt) a — n a /N, while the normal vector field is n a — 
g ab (dt) b /N, which involves the (s -I- 2)-metric, and as 
such, the lapse and shift. The arbitrareness of the lapse 
function and shift vector assures that one can still choose 
various foliations once n a is fixed by the chosen reference 
fluid. The time parameter t re f defined by the fluid as 
n a = (d/dt re f) a is different from the time t defined by 
the chosen foliation. Therefore while we associate t re f 
with the cosmological time, we still have the freedom 
of evolving the system with respect to any conveniently 
chosen time parameter t (like the conformal time). Re- 
stricting N = 1 and N a = leads to the identification 
t — t re j . 



The above conditions with Eqs. (|4Tj) and Eqs. (|43|) 
give 

Hf[g A ;V,p] = 



2n 2 M^ 



s-1 



H™[g A ;Va,Q a ] 
-2k 2 M^ 



s- 1 



V + (p + X)S(x) 



V a + Qa&(x) 



(49a) 



(49b) 



and 



SS 



M 



NM 



lab 



V~9 



s-1 



jpab 



+ ^V5[n Qb + (F-A).g ah ] ( 5(x), 

SS M ^ S 

SS M ^ S 

— = NV jj—v. 



(50a) 
(50b) 
(50c) 



The full brane-world geometrodynamics in the pres- 
ence of matter is then given by the equations (|2"9"]) for g A 
and the sum of the right hand sides of Eqs. ([30)1 and ([50]) 
for ir A : 

7r Qb = NS ab + NV ab — NMyfgC(L ab — Lg ab ) 

+^g~(D a D b N - g ab D c D c N - g ab D c ND c M) 
+V99 ab (d/d X - £m)(NC) + £ N ir ab 



-K 2 NM^/g 



_j,ab 



s-1 

+ [lP b + (P-A) g ab ]S( X )} , 
Pa = NVa-2^g-[L b a DbN + Da(N£)} 

;2 ; 



2s 



+ £ N Pa + K*Ny/<j -P Q , 

s — 1 

p = NS + NV - 2^fg(LL + D a D a N) 
s — 1 

The constraints (|42a[) and (|42b[) take the form 
-y/gM[R- L 2 + 3L ab L ab ) 



(51a) 
(51b) 

(51c) 



+ 



M y/g 



ab (d/d x 


— £m)L 


( „ ^ab 




1 Kab^ 












p a Pa + ~ 





pn 



+2k 2 M v / 3 



-V + (p + X)S(x) 



_s- 1 

0w H a [g A ,n A ;V a ,Q a ] = 

-2D b ir b a - (d/d X - £m)Pa+pD a M 
s 



, (52a) 



-2n 2 M^g 



s-1 



V a + QaS( X ) 



(52b) 
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The dynamical equations (|29|) and (|51j) and the con- 
straints (|52[) completely determine the time evolution of 
the geometry and the matter fields on the brane in brane- 
world scenarios. The LSDI matching condition follows 
from these equations, as we will show it in the next sec- 
tion. 

The 5-function type distributional sources in the evo- 
lution equations (|51a[) need some further interpretation. 
Such contributions also appear in the dynamics of K ab 
and /C, as derived in Paper I. These contributions indi- 
cate the singular behavior of ir ab . This is however, not 
surprising. The canonical equations and constraints are 
equivalent with the (s + 2)-dimensional Einstein equa- 
tions. If the sources of the latter are singular (in the 
present case across the brane), the Riemann (and Ein- 
stein) tensors are also singular, and certain singularities 
will be carried over in the canonical equations. Tradi- 
tionally (for example in the derivation of the effective 
Einstein equation [7[) the coefficients of the S- functions 
are interpreted as contributions present on the brane, but 
not in the bulk regions. 



V. REGULARIZATION ACROSS THE BRANE 

The brane contains (5-function type distributional 
sources, which in turn appear in both constraints (|52p 
and in the dynamical equation (|51ap . For these equations 
we apply the following regularization procedure. First we 
consider a domain of finite thickness (— \o, Xo) enclosing 
the brane at x = and we integrate the above-derived 
equations across its width. As a consequence, the Dirac 
distribution 5 (%) disappears. More precisely, according 
to [H, for any Tt (I) the relation 



dl6(l)H(l)=H(0) 



(53) 



holds. The integration across the brane is carried out 
here over a normal coordinate /, defined as d/dl = 1 
(see Q). Should we employ the coordinate x, defined as 
d/dx = Ml + M for integration across the brane (thus 
dl/dx = M), we obtain 



d X S (x)H(x) 



M 



(54) 



Secondly, the primitive function of the integral of any to- 
tal x-derivative term evaluated at the left and right do- 
main boundaries — xo and Xo give the so-called jump of 
the respective quantities. Thus, any quantity Q appear- 
ing as dQ / dx in the respective equation leads to its jump 
across the brane AC/, when the integration is carried out. 
Finally, if the value of xo is small, any other smooth func- 
tion of xo can be regarded as a constant, such that its 
integral will be proportional to the width 2xo ■ When we 
take the thin brane limit xo ~~ ¥ 0, this procedure drops all 
such terms, and what remains are only the terms origi- 
nally multiplying <5-functions and the jumps arising from 



the total x _ derivatives. If no such terms are present in 
any of the equations derived in the preceding sections, we 
obtain identities (of = type). Therefore non-trivial 
information arises only from the equations with total x~ 
derivatives and/or (^-functions. 

We can regain the junction condition for the embed- 
ding of the brane by integrating the constraints and those 
dynamical equations which contain derivatives with re- 
spect to Xi the coordinate running in the off-brane direc- 
tion. 

As Eqs. (|20ap and (|20b[) show, the vacuum constraints 
contain the x-derivatives of L ab and p a ■ The integration 
of the super-Hamiltonian constraint (|42a[) provides 



^g ab AL ab = -k 2 lim / d X M V? f T cd n c n d 

(55) 

By considering the brane world scenario, the integration 
of Eq. (|52a|) simply gives the trace of the junction con- 
dition (53) of Paper I: 



AL 



-K 2 (p + X) , 



(56) 



since /*; o d X f(x)S(x) = /(0)/M. 



When we integrate the supermomentum constraint 
(|42b|> . we obtain 



Ap a = 2k 2 lim / d X M^ (T cd g b c 7 



9c n d 9ab , (57) 



since the integral of the finite terms in the momentum 
constraint vanishes as xo —> 0. For the matter fields on 
the brane specified in Eq. (I49bj) . the integration of the 
momentum constraint (|52b[) gives 



Ap a = -2K 2 y/gQ a , 



(58) 



which is the vectorial projection (|B12[) of the Lanczos 
equation, rewritten in terms of momenta. 

By integrating the dynamical equation (|44a[) over Xj 
we obtain 

[AL afc - A(L - C)g ab ] 

d X M^T cd g a c g b d . (59) 

-Xo 



whereas the integration of the dynamical equation (|51a|) 
leads to the expression 

AL ab - A(L - C)g ab = -k 2 [H ab + (P - X)g ab ] . (60) 

After inserting Eq. (f56|) in the trace of this result, we get 



AC = k- 



2 (l-s)p-sP + X 



(61) 



which is the scalar projection (|B8|) of the Lanczos equa- 
tion. 

The substitution of Eq. (HTJ into Eq. gives 



AL ab = -k 2 n a6 + 



p + A 



-9ab 



(62) 
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Then we have obtained the tensorial projection (|B6[) of 
the Lanczos equation. This means the dynamical sys- 
tem ([5Tjl with the constraints (|52|) imply the usual LSDI 
junction conditions for the brane. 

By imposing Zi symmetry in the bulk across the brane 
(which implies AL a b = 2L a b, Ap a — 2p a , and AC = 2L), 
we can express the components of the extrinsic curvature 
associated with the brane normal in terms of the matter 
field variables: 



_ K I TT P + ^ 

L a b = I n a b H g ab 

2 \ s 

Pa = -^?\fgQa , 

_ „ 2 (l-s)p-sP + \ 
K 2s • 



(63a) 
(63b) 
(63c) 



As a simple application, we give the dynamical equa- 
tion of the heat flow. Eq. (|63b[) implies 



~2 i—f\ P a - be ■ 

K \fgQa = y 9 9bc - Pa 



(64) 



Here both gb c and p a are known as Eqs. (|29a[) and (|51h»|) . 
We obtain: 



Some of these equations contain 5-function type con- 
tributions, due to the singular source terms on the brane. 
These terms can be dropped, when we monitor gravita- 
tional dynamics in the bulk, and kept on the brane. 

The regularization of these equations across the brane 
yields the projections of the Lanczos equation, written in 
terms of canonical data. 

In the original covariant formulation of brane- world dy- 
namics Q the effective Einstein equation is obtained by 
expressing the terms quadratic in the extrinsic curvatures 
with matter variables. This is achieved by employing the 
Lanczos equation. In the present formalism, the role of 
these extrinsic curvatures are taken by L a b, £ a = fca and 
C, all functions of the canonical variables. The projec- 
tions of the Lanczos equation derived in this paper can 
also be employed to eliminate these geometrical variables 
in terms of matter variables in the canonical equations. 

As a simple application we have derived the equation 
of heat flow in terms of canonical data. 

The importance of the presented formalism relies in its 
possible application in the initial- value problem in brane- 
worlds and in the prospect of canonical quantization of 
brane-world gravity. 



K 2 y/gQa = Pa 



N 



-.p + D b N 



2V5 

-NV a + 2^[L\D b N + D a {NL)] 
2s 



-£-MPa - K Ny/g -P Q 

s — 1 



(65) 



This is the equation of heat flow expressed in terms of 
canonical data. 
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VI. CONCLUDING REMARKS 

We have derived the Hamiltonian dynamics of the 
(s + 2)-dimensional gravitation in terms of variables 
adapted to the existence of the preferred s-dimensional 
hypersurface. The canonical (gravi-tensorial, gravi- 
vectorial and gravi-scalar) metric variables g a b , M a arid 
M have canonically conjugated momenta 7r afc , p a and p, 
related to the extrinsic curvatures associated to the tem- 
poral normal n a . We have given the evolution equations 
for the canonical data and also the super-Hamiltonian 
and super-momentum constraints, all derived from an 
action principle. 



APPENDIX A: THE s + 1 + 1-DECOMPOSITION 
OF ENERGY-MOMENTUM AND THE 
EINSTEIN TENSORS 

In Appendix D of Paper I we have given the complete 
set of the decompositions of the Riemann-, Ricci- and 
Einstein tensors. Due to a typo, from the decomposition 
of the curvature scalar, Eq. (D3) of Paper I the term 
2 (K, 2 — C 2 ) was omitted, which cancels out the corre- 
sponding terms in the projections g^g b G c d, n a n b G a b and 
l a l b G ab - Therefore the formulae (D3)-(D4c) of Paper I 
correctly read: 



R 



R = R— 3K ab K ab + K 2 + 2 [{K + K)K + JC a JC a + g ab £ n K ab + £ n !C] 

+ZL ab L ab -L 2 + 2[{L- C)C - g ab £xL ab + £i£] 

-2 \N- 1 D a D a N + M~ 1 D a D a M + (NAiy 1 D a ND a M 



(Al) 
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and 



g c a g d b G cd = G ab -2K ac K c b + {K + K)K ab -2K a lC b + £ Yl K ab -N- 1 D b D a N 

\{ZK cd K cd - K 2 ) -{K + K)K - JCJC C - g cd £ Il K cd - £ n JC + N~ l D c D c N 



9ab 



2L ac L c b — (L — C) L ab - £iL ab - M~ x D b D a M 
^(3L cd L cd - L 2 ) + C(L -O- g cd £iL cd - £ X L - M~ X D C D C M g ab 
+ (NMy 1 D C ND C M , 

D c K c a — D a {K + /C) + K a L + CiJCa + M- l K c a D c M - M~ l KD a M , 
D c L c a — D a (L — £) + K a K + C n K a + N^L^D.N + N~ l LD a N , 

l - (R - K ab K ab + K 2 + 3L ab L ab - L 2 ) +KK- IC a JC a - g at ' £\L ab - M~ 1 D a D a M , 
~l(R + L ab L ab -L 2 - 3K ab K ab + K 2 ) - CL + K a K a - g ab £ n K ab + N~ l D a D a N 



9> d G cd 

g°J d Gcd 

n a n b G ab 

l a l b G a b 

n a l b G ab = D a K a - g ab £iK ab + K ab L ab + KL + M~ 1 IC a D a M 



(A2a) 
(A2b) 
(A2c) 

(A2d) 

(A2e) 
(A2f) 



The last three equations agree with Eqs. (D2d)-(D2f) of Paper I. From among them (|A2dp correctly yields the 
Hamiltonian constraint, Eq. (|13a|) . 



APPENDIX B: THE DERIVATION OF THE 
LANCZOS EQUATION 



We have demonstrated in Section V. that integration 
of the dynamical equations and of the constraints across 
the brane gives the projections of the Lanczos equation in 
the momentum phase space. In the velocity phase space 
these projections were obtained directly by projections of 
the full Lanczos equation, as Eqs. (53), (54), and (55) of 
Paper I (these are however valid only for s = 3). Equiv- 
alently, we can integrate the dynamical equations (67a), 
(67c) and the projection (D2e) of the Ricci tensor (i.e., 
the diffeomorphism constraint given in terms of velocity 
instead of momentum) of Paper I in order to obtain in 
full generality (for generic s) the projections of the Lanc- 
zos equation in the velocity phase-space. With this we 
fully establish the commutativity of the variation princi- 
ple and the geometrical decomposition of the quantities, 
to be discussed in more detail in Appendix C. 

In order to carry out this program, we apply the regu- 
larization procedure described in Section V. for the evo- 
lution equations of K ab and /C, derived in Paper I, which 
contain the following projections of R a b- 



rewritten conveniently as 



9 K 
di Kah 



k 2 N 



(p + X) 



g a b + n a . 



d_ 

Of 



9 ( N T \ A 

+ \M Lab ) +Aab ' 
K = £ 2 A>-^ + (s + 1)A] 



d_ 

dx 



Hx) 

(B3a) 

Hx) 

A , (B3b) 



where in A ab and A we have collected only finite terms: 
2. 



N 



-K 2 Ag ab - R ab + L ab (£-£)- 2L ac L c b 



-K ab (K + JC) + 2K ac K c b + 2K a lC b 
~ (D b D a M ~ M c D c L ab - 2L c{a D b) M c 
-D b D a N + N c D c K ab + 2K c(a D b) N c 

(^) , (B4a) 



glgtRcd = k 2 

l a l b R ab = k 2 



2A+(p + X)S( X ) 



g a b + n Qb <5 (x) 



2A+[p- sp+(s+l)X]8(x) 



(PI) 
(B2) 



The first and third of the Eqs. (67) of Paper I can be 



A = Nl-K 2 A-L ab L ab + £ 2 -2K a K a -K{K + K) 



1 „x D a D a M 

■m md ^ l -^ + ^aT- 



D a M 
M 



D a N + N a D a K 



(B4b) 
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The integration across a finite coordinate distance con- 
taining the brane (for example from — xo to Xo) °f Eqs. 
(|B3|) and the subsequent limit xo — ► gives 



:2 N 

" M 

,N_ 

M 



(p + X) 



-dab 



n„ 



A I Tl L ab 



p — sp + (s + 1) A 



N 
M 



(L-C) 



(B5a) 



(B5b) 



(In deriving Eqs. (|B5[) we have employed that time- 
derivatives and integration over x commute.) With the 
remark that neither M, nor N are discontinuous across 
the brane, but the extrinsic curvature L ab is (as it de- 
pends on the brane embedding into the bulk on each 
side), from Eq. (|B5aP we obtain the jump of L ab as 



AL ab = 



(P + A) 



g a b + n at 



(B6) 



This is the tensorial projection of the Lanczos equation, 
Eq. (53) of Paper I, valid for generic s. The trace of the 
right hand side of this equation is 



AL = ~k 2 (p + X) . 
With this, Eq. (|B5bj) implies 



(B7) 



(B8) 



which is the scalar projection of the Lanczos equation, 
Eq. (55) of Paper I, valid for generic s. 

The vectorial projection of the Lanczos equation does 
not emerge in a similar way. Indeed, the remaining evo- 
lution equation for JC a [the second Eq. (67) of Paper I], 
is 



dt Ka 



N [-D b L ab + D a {L-C)- KfC a ] + N b D b K, a 
- (L b a + C6 b a ) D b N + !C b D a N b . (B9) 



As all terms are regular, integration over the range 
( — Xo, Xo) an d the limit xo - * would give nothing but 
the identity = 0. 

The vectorial projection of the Lanczos equation can 
be instead obtained from Eq. (D2e) of Paper I, which 
contains the projection g c a n d R cd , given by the bulk Ein- 
stein equation: 



g c a n d R cd = k 2 



g c a n d U cd ~ Q a 6 (x) 



Written in terms of x-derivatives, Eq. (D2e) of Paper I 
becomes 

= K 2 Q a S (x) + -jjj -^ICa + A a (BIO) 



K 2 g c a n d U cd + D c K r a -D a (K + lC)+ K. a L 



(KlDiM 



KD a M + K c D c M a - M c D c JC a 



where A a represents the collection of finite terms. Inte- 
gration over an infinitesimal range (— Xo, Xo) and taking 
the limit xo — > gives 



AIC a = ~K 2 Q a 



(B12) 



This is the vectorial projection of the Lanczos equation, 
Eq. (54) of Paper I. 



APPENDIX C: DYNAMICS REGAINED FROM 
THE (s + 1) + 1 ADM DECOMPOSITION 

In this Appendix we regain the full (s + 1 + ^-break- 
up of the (s + 2)-dimensional equations by a step-by- 
step method, by further splitting of the ((s + 1) + 1)- 
dimensional ADM decomposition of the bulk. First we 
foliate the (s + 2)-dimensional space-time B with the 
(s + l)-spaces St- Then we derive the equations of mo- 
tion together with the constraints from the corresponding 
variation principle. Finally we embed the s-spaces E* x in 
the hypersurfaces S t and perform a second split-up with 
respect to the extra dimension of the geometrodynamics. 

We begin by briefly presenting the standard (s + 1) + 1 
ADM decomposition of the vacuum Lagrangian. This 
yields: 

C G = N^{R-k 2 + K ab k ab )-2V a [N^{a a -kn a )] , 

(CI) 

with the dynamical variables 'gab and K a b, representing 
the first and second fundamental forms of St, respec- 
tively. Here R is the intrinsic curvature scalar of St and 
K = g ab K a b is the trace of K a i,. This equation can be 
rewritten in the form 

+N a H G - 2^D a a N + N b k ab )} (C2) 

containing the pure spatial projections £tK ab = 
g c a g d b£tK cd and D a K bc = g d a Tbg f c ^dK ef of the 
(s + 2)-dimensional Lie-derivative £tK cd and covariant 
derivative V d K e f, respectively. We can also identify the 
cofactors of the Lagrange multipliers N and N a as the 
super-Hamiltonian and super-momentum constraints 

Hi = -2^n a n b G ab = Vd(R + K 2 -k ab k ab yV3ei) 
H G = -2VW b a n c G bc ^2VdD b (k ab - k<f b ) ■ (C3b) 
By introducing the (s + l)-dimcnsional momentum 

dL G 



d(£tdab) 



g(K ab -r"K), (C4) 



(BH) th e vacuum action S G — J d s + 2 x^/ —gL G can be cast 
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into the form 

S G [g ab ,n ab ;N a ,N] 



l dt ) 


f d s ^ 




N a H G ) 


-\ 


[r 


^ 2 x^/-^V a {a a - Kn a ) 




- 2 J 


[ dt 


f d s+1 xD a (N b w ab ) , 


(C5) 



and the constraints transform to 
'H ( l(x;g ab ,n ab ] = -y/§R 



1 V 



(C6a) 
(C6b) 



H a (x;g ab ,n ab ] = -2D b 9 ab 



for any point x £ St- Extremizing the action (|C5|) with 
respect to the canonical variables 'gab and n ab yields the 
field equations 



£tg a .b 

X t 7T 



2Ng 



-1/2 



TTab 7??ab ) + ^Nffafc , (C7) 

n 



-2Ng 
—Ng 



-1/2 I -a ~bc 

1 1 7T r lT 



-1/2 



1 



ab 



-Ng 1 / 2 (^R ab - ]^Rg ab 
-g 1/2 {D a D b N -g ab D c D c N) , 



(C8) 



whereas extremizing the action with respect to the 
Lagrange multipliers N and N a yields the super- 
Hamiltonian and the super-momentum constraints of 
vacuum gravity: li. G — and H a = 0. 

The dynamical equations (|C7|) and (|C8|) can be equally 
interpreted as the time evolution of the canonical vari- 
ables, generated by the smeared constraints H G [N] = 
H G [N}+ H G [N a }: 

£tg ab = {g ab ,H G [N]} , (C9a) 
£ t 9 ab = {9 ab ,H G [N}} , (C9b) 

where the Poisson bracket is defined by 

{f(x),h(#)} = I d^x» * f 

J s t og ab {x") 6ir ab (x") 

_ f Sf(x) 5h(x') 

for any function f(x;g ab ,TT ab ] and g(x; g ab , ir ab ]. The 



Poisson brackets of the constraints give the Dirac algebra 

{H G (x),H G (x')} = g ab (m G (x)6, b (x,x') 

-{x<r+x') , (Clla) 

{H G (x),H G (x>)} = H G ±(x)S, a (x,x') 

+H G _Jx)S(x,x') , (Cllb) 

{H G (x),H G (x')} = H G (x)5, a (x,x') 

- (ax <-> bx') . (Cllc) 



The dynamical equations (|C7[) and (|C8[) together with 
the constraints (|C6[) provide the full geometrodynamics 
of the vacuum gravity. Any matter field can be equally 
((s + 1) + l)-dimensionally decomposed and coupled to 
gravity. 

The next step necessary in order to recover the 
(s + 1 + l)-dimcnsional decomposition of geometrody- 
namics is a further split. We can interpret the momenta 
7r ab , p a and p introduced in the main text as the projec- 
tions of Ti ab : 

-ab = n ab + Ml (a p b) + M-tl b p . (C12) 



By inserting Eq. (|C12[) into the action (|C5|) . we re 
gain the fully decomposed action Q18p with the con- 
straints (|20|) . Similarly, the dynamical equations (|C7| 
and (|C8p can be split by applying the formula (|C12[) 
The only non-trivial step is the decomposition of the Lie- 
derivatives. For this we employ: 



(£ t — £^)Tr ab 



M,„, h { d 
2 

N fs-l 



Mp 2 - np ) l a l b . (C13) 



2s^g V 2 



With this we regain the fully decomposed dynamical 
equations (j2l?a) - (PDc]) . From Eqs. (|C12|) and jC6| we 
can also obtain the constraints (|2U|) . by using the twice 
contracted Gauss equation for the hypersurface Y> tx of 
S t - 

R = R + L 2 - L ab L ab + 2D a (l c DJ a - Ll a ) . (C14) 

We have thus shown that the variation principle and 
the (s + 1 + l)-decomposition commute, i.e. no matter 
which order we apply the decomposition of the bulk and 
the extremization of the action, the result is the same. 
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